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We examine how corrections to S-state energy levels, EnS, in hydrogenic atoms due to the finite 
proton size are affected by moments of the proton charge distribution. The corrections to EnS are 
computed moment by moment. The results demonstrate that the next-to-leading order term in 
the expansion is of order r^jas times the size of the leading order (rp) term. Our analysis thus 
dispels any concern that the larger relative size of this term for muonic hydrogen versus electronic 
hydrogen might account for the current discrepancy of proton radius measurements extracted from 
the two systems. Furthermore, the next-to-leading order term in powers of r^jas that we derive 
from a dipole proton form factor is proportional to (r^), rather than (r^) as would be expected 
from the scalar nature of the form factor. The dependence of the finite-size correction on (r^) and 
higher odd-power moments is shown to be a general result for any spherically symmetric proton 
charge distribution. A method for computing the moment expansion of the finite-size correction to 
arbitrary order is introduced and the results are tabulated for principal quantum numbers up to 
n = 7. 


I. INTRODUCTION 


The proton radius puzzle concerns the incompatibly 
different measurements of the proton radius 
resulting from normal hydrogen (ep) spectroscopy and 
spectroscopy of muonic hydrogen {pp). The 2010 recom¬ 
mended CODATA value based on hydrogen spectroscopy 
and electron-proton scattering is Vp = 0.8775(51) fm, 
[I| whereas the CREMA collaboration found rp = 
0.84087(39) fm from precision measurements of the Lamb 
shift in pp. ii In atomic spectroscopy, the size of (r^) 
is calculated from its shift of the A-state energy levels, 
which is given to leading order in powers of the proton 
radius by 


AEr^s = ^ {tD = ( 1 ) 

In scattering experiments, Vp is obtained from the Sachs 
electric form factor Ge{Q'^) according to 



dGsiQ^) 

dm 


Q2=0 


( 2 ) 


Reviews of the proton radius puzzle by R. Pohl et al. 
0 and C. Carlson @ discuss in detail the importance of 
this problem and include various theoretical explanations 
that have yet to gain acceptance. 

Among the explanations of any discrepancy between 
theory and experiment, there must always be included 
the possibility that the theory was not correctly applied, 
and in this article we investigate one particular possibil¬ 
ity. We note that Eq. ([21) arises from using the Sachs 
electric form factor Ge{Q^) in a first-order expansion in 
terms of As pointed out in Q, the proton ra¬ 

dius puzzle makes it natural to consider the effects of the 
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and higher-order terms in this expansion. Retaining 
such terms in the form-factor and proceeding as with the 
first-order calculation leads to the proposal of a lepton- 
dependent parametrization of the leading order finite-size 
energy correction to the nS levels of hydrogenic atoms, 

AEns,e = ^IV’nS,£(0)pr2_^, (3) 

where £ € {e,/i}. Since the muonic Lamb shift measured 
by the CREMA collaboration concerns a 2S-2P tran- 
sitioin we demonstrate that the proposed parametriza¬ 
tion [6| of Tp ^ leads to a divergence in the energy correc¬ 
tion of Eq. ©■ We go on to examine the argument of 
@ based on continuing the expansion of Eq. in pow¬ 
ers of and confirm that the next term in the series is 
indeed more significant for pp than for ep. However, as 
will be shown, the effect remains far too small to account 
for a 4% discrepancy in rp. Our analysis parallels the 
work of 0, Hi , which reached the same conclusions using 
numerical methods, and we generalize those results. 


II. DIVERGENCE OE THE PROPOSED 
CORRECTIONS 

The proposed lepton-dependent radius extracted from 
nS energy levels which appears in Eq. m arises from 
expanding in powers of q^ 0], an idea that is well- 
motivated. Keeping terms up to q'’’ and rearranging the 
result so as to match the factorization in Eq. m would 
seem to suggest the reparametrization, 

= ’’p - {qL^)n,efp, (4) 

(q^)n.f = |^/’nS,^(0)r^y j 

( 5 ) 

with Pp a parameter determined by the second derivative 
of Ge. An explicit expression for will not be needed. 
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Let us now consider as an example the 2S wavefunction, 
which is well-known to be 

where 

ai = (7) 

GlTTZred,^ 

denotes the Bohr radius of either ep or pp, and we work in 
center-of-mass coordinates. Inserting this wavefunction 
into Eq. yields 


correction to the nS level can be written as 

SEr^s = I d\\^^s{r)f I ( 10 ) 

where (5E(q) is the Fourier transform of Eq. Q, 

<5E(q) = ^ [1 - E(q)] , (11) 

with F(q) being the NR form factor. To proceed with 
analyzing the higher-order terms in 5E„s we first approx¬ 
imate the form factor with what is traditionally known 
as the dipole form factor: 



Evaluating the q-space integral first leads to an im¬ 
mediate divergence. Equivalently, using the fact that 
|q| 2 giq.r _ _y 2 giq.r integrating by parts, one finds 
that this divergence corresponds to trying to evaluate 
V^|^/’ 2 S,f(r)P at the origin, which can be seen to diverge 
using the Schrodinger equation. These considerations, 
which hold for any S'-state, show the proposed finite-size 
correction is divergent, so the idea given by @ cannot 
bring the different measurements of Vp into agreement. 
Nevertheless, it is worthwhile to pursue the idea that 
higher q^" terms might be important for the finite-size 
correction in muonic hydrogen. This is done in the fol¬ 
lowing by first evaluating the relevant matrix element as 
a function of and then expanding the result in terms 
of the ratio of the size of the proton to the Bohr radius. 
This calculation will therefore show that using the proper 
order of integration and series expansion is essential to 
getting the correct result. 


III. EXPANSION OF THE S'-STATE 
FINITE-SIZE CORRECTION 


We now examine what the true effect of higher-order 
moments is on the size of the nS-level corrections derived 
from nonrelativistic (NR) perturbation theory. 


A. Calculation from an approximate NR proton 
form factor 


The perturbation to the Hamiltonian is given in coor¬ 
dinate space by 


SV{r) 



>(r') -(5(r') ~ 


(9) 


where p describes the electric charge distribution of the 
proton (normalized to unity). The corresponding energy 


f7ii:(q) — Edipoie(q) 



( 12 ) 


where ~ 0.71 GeV^. The name dipole arises from the 
fact that Ge has a second-order pole at |qp = —A^. Sub¬ 
stituting the dipole form factor into Eq. m and evalu¬ 
ating the q-space integral of Eq. (fTH leads to 


n / Ar -I- 2\ 

SE^S =aj^ drr^ |i?„o(r)|" j 

In particular, one finds that 


5E2S = olK 


4(a^A) ^ -|- 3(a£A) ® -|- (a^A) ® 
4(1 -I- (a^A)-!)^ 


(14) 


Since (a^A) ^ 1 for both ep and pp^ this quantity 

can be used as an expansion parameter: 


dE2S = a A 


(ofA)”^ - 5(a^A)“^ -I- —(a^A)”'" 


77, 


-y(a^A)-'> + 80(a^A)-"- 


(15) 


It will be economic to define ^nS {p, as the ratio of the 
second term to the first term in an expansion of the form 
in Eq. m for SEnS, resulting from the distribution p 
with satellite £. In order to account for a 4% reduction 
in the apparent rp (an 8% reduction in r^), we require 
^ 2 s(Pi p) ~ 8%. Using the known values of ai and A, we 
calculate 


^2s{p,£) 


-4.1x10-3, £ = p, 

-2.2x10-3, £ = e, 


(16) 


for a dipole form factor. The relative size ^ 2 s{PtP) = 
—4.1 X 10-3 therefore shows that neglecting higher-order 
moments in Eq. © is not to blame for the proton ra¬ 
dius discrepancy [despite the fact that ^nS {p, £) is about 
200 times greater for pp versus ep]. More concretely, as 
discussed in 0 , any missing terms from the energy differ¬ 
ence associated with the muonic Lamb shift would need 
to amount to 0.31 meV in order to bring their value of rp 
into agreement with the previously accepted CODATA 
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value, and the second term of Eq. dm amounts to only 
a 0.02 meV difference. The conclusions we have reached 
so far, namely the inadequacy of higher-order terms in 
the finite-size correction as a solution to the proton ra¬ 
dius puzzle, will next be substantiated without assuming 
a specific shape of the proton form factor. 


B. General dependence of E2S on moments of the 
proton charge distribution 

In the preceding section, the assumption of a dipole 
form factor lead to the expansion in Eq. Uni) for 5E2S- 
However, Eq. US) can be cast in a more revealing form 
by trading the powers of A for moments of the proton 
charge distribution, 

(r;)=|p(r)r"dV (17) 

The dipole form factor corresponds to a proton distribu¬ 
tion function, 


p(r) = A^c-^VStt, 

(18) 

from which one finds that 


A- = 2(r;)/(n + 2)!. 

(19) 

Thus, the expansion in Eq. (11511 reads 


.. , YD 

12a| 120^ 160a| 

iHrl) (rp 

720a^ 252aJ 

(20) 

At this point we call attention to the peculiar result 
that all moments of Vp higher than the first appear in 
the finite-size correction Eq. (1^01). It is straightforward 
to show that the form factor of a spherically symmetric 
density is given by 


Jiff ( — 

-P'(q)lp(r)=p(|r|) = X! (2m + 1)! 

i.e., the form factor can be written purely in terms of even 
moments of the coordinate-space density. In light of Eq.. 
m and (ED, the appearance of odd-power moments in 
Eq. (I^Dl) is surprising. We believe that the form of the 
series in Eq. (I^Tl) may have been the motivation for the 
lepton-dependent parametrization of Eq. O- We note 
that the term proportional to (r^) in Eq. (I2L)H i s not at 
all related to the third Zemach moment [2, [l^, which 
produces a much larger change to the energy level. 

Although we obtained the expansion of Eq. (1^ by as¬ 
suming a dipole form factor, this is in fact the general 
result for any spherically symmetric proton charge den¬ 
sity. Other models for the density which can be verified 


to yield the same expansion as Eq. (12011 include 


PG(r) = 


27r(r2) 


p/ J 


3/2 


exp - 


3r^ 

W^) 


PY(r) 


3 / -V6r \ 

27r(r2)r 


( 22 ) 

(23) 


although it should be noted that the values of the mo¬ 
ments themselves depend, of course, on which density is 
used. (Densities of the above forms were also used in the 
numerical work of Q.) 

To justify our result in general, we claim that the 
spherically symmetric density appearing in Eq. ([HD— 
whatever its particular form may be—can be represented 
formally by the series 


P(i’) 


1 

47rp2 


00 


E 


j! 




(24) 


where S(r)dr = 1 and S^^\r) = {-^yS{r). This al¬ 
lows us to express the finite-size correction to the nS level 
as 


SEnS = -^E 

i=i 


dr'r' \Rnoir')y / dr 


i-YH) 

j'-r> 


SU)(r) 


(25) 

where r> = max(r, r'). 

It is not immediately obvious how one should go about 
evaluating the terms in Eq. (1^51) since S{r) lives on the 
non-negative real axis rather than the usual entire real 
line; typical methods such as integration by parts lead to 
ambiguous expressions involving, e.g., <5(0) or (5'(0). One 
solution in such circumstances is to replace the Dirac 
delta with a representative test distribution /(r, e) which 
becomes S(r) in the limit e —> 0. However, many of the 
usual test functions diverge for terms j > 2 when the 
e —>■ 0 limit is taken. 

We find that the following sequence of test functions 
can reliably be used to evaluate any term of the series in 
Eq. (1^ : 


fjir, e) 


e /r\t+2 

{j + 2)!e Ve/ 


(26) 


These test distributions are unit-normalized on the posi¬ 
tive real axis for any finite e > 0 and the jth test function 
is used to evaluate the jth term of Eq. (1^ . In contrast 
to the usual prototypes used for the entire real line, our 
test functions peak on the positive real axis and only 
reach the origin as e —?► 0. 

Evaluating the terms of Eq. using the test func¬ 
tions described above and taking e —>■ 0 at the end yields 
the anticipated expansion for SE 2 S in moments. 


_ , jrl) (4) , Urt) 

12a^ 12a^ 160af 

11 ^^) K) 

720a^ 252aJ 


( 27 ) 
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Table I: Numerical coefficients Cn ^ of the moment expansion for SEnS (cf. Eq. (1271) 1. The coefficients are calculated according 
to the method described in section HVl The j=l column of zeroes is just the restatement that the finite-size correction to E^s 
has no dependence on (rp). The j = 2,3 columns suggest that Cn' = —Ck for all n. 


Cn^ i=l j=2 j=3 j=4 j=5 j=6 j=7 j=8 j=9 j=W 


n= 

T 

0 

2 

-2 

2 

-8 

4 

-2 

2 

-16 

4 

3 

3 

5 

45 

63 

105 

405 

14175 

17325 

n= 

-2 

0 

1 

-1 

7 

-11 

1 

-11 

29 

-37 

23 

12 

12 

160 

720 

252 

13440 

207360 

1814400 

8870400 

n= 

Q 

0 

2 

-2 

46 

-136 

404 

-122 

206 

-7888 

3004 


81 

81 

3645 

32805 

413343 

688905 

7971615 

2511058725 

9207215325 

n= 

4 

0 

1 

-1 

27 

-13 

11 

-113 

1891 

-3679 

421 


96 

96 

5120 

7680 

28672 

1720320 

212336640 

3715891200 

4541644800 

n= 

:5 

0 

2 

-2 

42 

-8 

1556 

-1294 

3254 

-60272 

162556 

375 

375 

15625 

9375 

8203125 

41015625 

791015625 

138427734375 

4229736328125 

n= 

6 

0 

1 

-1 

181 

-257 

89 

-1553 

9161 

-74143 

92579 

324 

324 

116640 

524880 

826686 

88179840 

4081466880 

321415516800 

4714094246400 

n= 

■7 

0 

2 

-2 

82 

-232 

3476 

-134 

155894 

-550288 

391516 

1029 

1029 

84035 

756315 

51883209 

12353145 

114402475845 

4004086654575 

34257185822475 


Thus, taking the ratio of the first two terms, we have the 
general result 


Us{pJ) = -p^-^ (28) 

ae{rp) an 

as long as p has a range much smaller than the Bohr 
radius an- Note that the generalization from 2S to nS 
is inferred from Table [H We are now able to claim— 
without any assumptions about the specific shape of the 
proton charge density or, equivalently, the form factor— 
that the inclusion of higher-order moments in Eq. m is 
of negligible impact on the value of (r^) inferred from 
measurements of the 2S-2P Lamb shift. 


IV. MOMENT DEPENDENCE OF THE 
FINITE-SIZE CORRECTION FOR n / 2 


Writing the finite-size correction Eq. ((Mil in the form 


OO 



. 0 ) 


(ri) 


(29) 


we can summarize the dependence of EnS on moments 
of the proton charge distribution via the dimensionless 
coefficients Cn ^. In Table (H we list the coefficients for 
n = 1,...,7 and j = 1,...,10 as calculated by the 
method introduced in the previous section. The tabu¬ 
lated coefficients suggest that = —Cn^ for all n, so 
that Eq. (IMl) holds for any S'-state. 


ton radius measurements obtained from electronic versus 
muonic hydrogen. Using a typical approximation to the 
proton form factor as well as a more general coordinate- 
space analysis, we have determined the size of the effects 
of higher-order moments of the proton charge distribu¬ 
tion. The expansion of Eq. (1291) combined with the coef¬ 
ficients in Table H] could also be used for nuclear targets 
for which the sizes of the moments are more significant, 
as well as for transitions involving higher S'-states in light 
atoms. Our calculations Eq. (fMl) and Eq. (IMl) show that, 
while it is true that the next term in a series expansion 
for the finite-size correction in powers of rpjan is more 
important with muonic hydrogen, the correction cannot 
be the source of a 4% discrepancy in rp. Furthermore, 
this term is not 0{{rp)), but is ©((r^)). The finite-size 
correction to Ens depends quite generally on all moments 
of the proton charge distribution of at least second-order. 
This unexpected result has been verified using several dif¬ 
ferent approaches. We close with the reaffirmation that 
truncating the finite-size energy correction Eq. ([T]) at (r^ 
can make a difference of no greater than a few parts in 
10^ for p,p and that the proton radius puzzle remains 
open in this regard. 
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